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Introduction 

Earth  monitoring  systems  of  the  future  may  include  large  numbers  of  inexpensive  small  satellites,  tasked  in  a 
coordinated  fashion  to  observe  both  long  term  and  transient  targets.  For  best  performance,  a  tool  which  helps 
operators  optimally  assign  targets  to  satellites  will  be  required.  We  present  the  design  of  algorithms  developed 
for  real-time  optimized  autonomous  planning  of  large  numbers  of  small  single-sensor  Earth  observation 
satellites.  The  algorithms  will  reduce  requirements  on  the  human  operators  of  such  a  system  of  satellites, 
ensure  good  utilization  of  system  resources,  and  provide  the  capability  to  dynamically  respond  to  temporal 
terrestrial  phenomena.  Our  initial  real-time  system  model  consists  of  approximately  100  satellites  and  large 
number  of  points  of  interest  on  Earth  (e.g.,  hurricanes,  volcanoes,  and  forest  fires)  with  the  objective  to 
maximize  the  total  science  value  of  observations  over  time.  Several  options  for  calculating  the  science  value  of 
observations  include  the  following:  1)  total  observation  time,  2)  number  of  observations,  and  the  3)  quality  (a 
function  of  e.g.,  sensor  type,  range,  slant  angle)  of  the  observations.  An  integrated  approach  using  integer 
programming,  optimization  and  astrodynamics  is  used  to  calculate  optimized  observation  and  sensor  tasking 
plans. 

Problem  Motivation 

“...Thus  far,  we  are  only  experimenting  with  long  term  weather,  climate,  and  natural  hazard  prediction.  The 
quest  for  a  true  predictive  capability  for  Earth  system  changes  requires  a  flexible  and  progressive  space  system 
architecture  that  is  responsive  to  our  needs  based  on  our  current  understanding  of  the  system  as  well  as 
accommodating  emerging  needs  in  the  coming  decades.  We  need  to  design  and  establish  a  smart,  autonomous 
and  flexible  constellation  [of]  Earth  observing  satellites  which  can  be  reconfigured  based  on  the  contemporary 
scientific  problems  at  hand.  Such  a  constellation  would  exploit  a  combination  of  active  and  passive  sensing 
sensors  in  ways  that  we  can  perhaps  imagine  today . ” 

This  is  a  quote  from  the  remarks  of  former  NASA  Administrator  Daniel  S.  Goldin  “The  Frontier  of 
Possibilities”  presented  at  the  International  Astronautical  Federation  on  October  3,  2000.  It  clearly  describes 
the  underlying  rationale  for  the  multi-year  project  in  which  we  are  developing  algorithms  for  resource 
allocation  via  autonomous  reconfiguration  of  satellite  webs  consisting  of  heterogeneous  Earth  observation 
sensor  platforms. 

The  development  of  these  algorithms  and  a  simulation  testbed  in  which  they  reside,  the  Earth  Phenomena 
Observing  System  (EPOS),  will  reduce  requirements  on  the  human  operators  of  satellites,  improve  system 
resource  utilization,  and  provide  the  capability  to  dynamically  respond  to  temporal  terrestrial  phenomena. 
Examples  of  triggering  events  are  localized  transient  phenomena  that  have  a  significant  impact  on  human  life 
such  as  volcanic  eruptions,  weather  (hurricanes,  tornadoes,  etc.),  algae  plumes,  large  ocean  vortices,  ice  shelf 
break-up,  seismic  activities,  oil  spills,  magnetic  anomalies,  and  search  and  rescue. 


Paper  presented  at  the  RTO  AVT  Course  on  “Intelligent  Systems  for  Aeronautics” , 
held  in  Rhode-Saint-Genese,  Belgium,  13-17  May  2002,  and  published  in  RTO-EN-022. 
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Vision 

In  the  year  2020,  automated  mission  management  will  be  required  to  help  NASA  achieve  their  Sensorweb 
vision.  Draper’s  vision  is  to  provide  functionality  for  that  mission  manager,  EPOS  2020.  Key  technology 
required  for  the  full  EPOS  2020  is  being  developed  under  our  current  2  to  3  year  effort. 

In  2020,  many  Earth  observing  systems  will  be  in  place,  each  with  its  own  ground  facilities,  communication 
and  data  protocols,  satellite  characteristics  and  observation  schedules.  There  will  be  hundreds  of  Earth 
observing  satellites  and  relatively  few  will  be  identical.  In  most  systems,  satellite  maneuvering  will  be 
reserved  for  station-keeping  only.  From  EPOS  perspective  these  are  "coasting  satellites"  -  EPOS  observation 
schedules  must  utilize  the  satellite’s  existing  orbit.  Some  of  the  systems  have  satellites  that  can  be  dynamically 
tasked  through  in-plane  maneuvering.  For  EPOS  purposes  these  are  referred  to  as  "maneuvering  satellites." 
These  satellites  will  be  easily  upgraded  and  refueled  using  technology  currently  being  developed  [5].  We 
hypothesize  there  will  be  a  coordinating  ground  station  that  has:  1)  communication  with  these  Earth  observing 
systems,  2)  a  given  set  of  data/information  interchange  protocols  that  provides  the  mechanism  to  exchange 
data  with  these  varied  systems,  3)  visibility  into  whatever  parameters  are  needed  to  model  these  systems,  and 
4)  the  capability  to  influence  the  observation  schedules  of  these  systems.  This  is  illustrated  in  Figure  1.  When 
a  target  becomes  known,  the  coordinating  ground  station  calculates  which  of  the  satellites  controlled  by  the 
Earth  observing  systems  can  provide  suitable  information  at  the  appropriate  times  and  uses  that  knowledge  to 
send  appropriate  viewing  requirements  to  the  individual  system’s  ground  stations.  The  individual  ground 
stations  look  more  closely  at  the  availability  and  appropriateness  of  the  satellites  within  its  control,  especially 
dealing  with  the  priority  of  all  requests  for  satellite  assets.  An  inability  to  satisfy  the  viewing  requirements  is 
communicated  back  to  the  coordinating  ground  station. 


Figure  1:  Year  2020  Operation  Concept 


This  concept  of  operations  includes  planning  and  control  of  constellations  and  collaborative  groups. 
Constellations  are  defined  as  a  system  employing  two  or  more  spacecraft  whose  orbits,  operations,  and 
observations  are  coordinated  to  provide  global  coverage  or  to  improve  temporal  resolution  from  an  altitude 
below  GEO  [4],  Constellations  operate  within  constraints  (e.g.,  keeping  relative  position  in  an  orbital  plane) 
that  a  planning  and  control  system  must  recognize.  Collaborative  groups  are  temporary  sets  of 
satellites/virtual  platforms  that  have  commonality  in  their  objectives;  satellites  can  simultaneously  belong  to 
two  collaborative  groups.  These  groups  “exist”  (i.e.,  are  formed)  in  the  coordinating  ground  facility  for  the 
time  period  during  which  they  are  executing  a  plan.  In  our  concept  of  operations,  a  satellite  within  a 
constellation  can  also  be  a  member  of  a  collaborative  group  as  long  as  the  constellation  constraints  are  not 
violated. 
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Approach 

The  implementation  of  EPOS  combines  Draper’s  autonomous  system’s  planning  framework  [3],  which  is  used 
to  decompose  the  problem  into  tractable  subproblems  or  levels.  Each  level  is  then  solved  with  a  combination 
of  combinatorial  optimization  combined  with  astrodynamic  modeling. 


Planning  Framework 

The  basic  building  block  of  Draper’s  planning  and  control  architecture  for  autonomous  systems  [3]  is  shown  in 
Figure  2.  It  is  an  extension  of  the  sense-think-act  paradigm  of  intelligence,  and  is  similar  to  the  military’s 
Observe-Orient-Decide-Act  loop  [2] .  The  key  elements  are  modules  for  situation  assessment,  plan  generation, 
plan  implementation,  and  coordination 


■  Supervision  from  Higher  Levels 

-  objectives 

-  constraints 

’  Information  from  Cooperative  Agents 
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Figure  2:  Planning  and  Control  Architecture  Modules 

The  planning  problem  addressed  in  our  effort  is  complex  enough  to  warrant  hierarchically  decomposing  it  in 
order  to  make  it  tractable.  Replanning  takes  place  at  the  lowest  level  possible,  without  disturbing  other  plans 
unless  additional  resources  are  needed.  Hierarchical  decomposition  is  appropriate  for  applications  in  which 
there  is  significant  stochasticity.  In  such  systems,  it  is  not  practical  to  make  detailed  plans  too  far  into  the 
future,  since  the  state  of  the  world  (e.g.,  new  targets  of  interest)  and  the  state  of  the  system  (e.g.,  a  satellite 
sensor  fault)  can  change. 

In  order  to  make  the  optimization  problem  tractable,  hierarchical  decomposition  is  used  both  temporally  and 
functionally.  The  decomposition  is  characterized  by  higher  levels  that  create  plans  with  the  greatest  temporal 
scope  (longest  planning  horizon)  but  with  the  least  detail.  At  lower  levels,  the  planning  horizon  becomes 
shorter  (nearer  term),  but  the  level  of  detail  of  planned  activities  increases.  The  less  detailed  plans  at  the 
higher  levels  coordinate  or  guide  the  generation  of  solutions  generated  at  the  lower  levels.  Indeed,  planning 
actions  over  extended  periods  of  time  at  a  high  level  of  detail  is  typically  both  futile  and  impractical.  Futile 
because  detailed  actions  planned  on  the  basis  of  a  specific  prediction  of  the  future  may  become  obsolete  well 
before  they  are  to  be  executed  due  to  an  inability  to  accurately  predict  the  future.  Impractical  because  the 
computational  resources  required  to  develop  detailed  plans  over  extended  periods  of  time  may  be  prohibitive 
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either  in  cost  or  availability  or  both.  The  relationship  between  the  levels  of  the  hierarchy  and  the  planning 
horizon  and  level  of  plan  detail  is  shown  in  Figure  3. 
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Figure  3:  Temporal  Decomposition 

One  way  of  grouping  the  types  of  decisions  that  need  to  be  made  for  EPOS  is  shown  in  Figure  4.  Three 
decision  tiers  are  present:  1)  System,  the  top  tier  focuses  on  decisions  that  impact  the  entire  EPOS  (e.g.,  which 
targets  to  collect  data  for);  2)  Collaborative  Group,  the  middle  tier  addresses  issues  that  relate  to  a  group  of 
multiple  satellites  being  used  to  collect  data  on  a  target  (e.g.,  which  satellites  make  up  this  group);  and  3) 
Satellites,  the  lowest  tier  addresses  decisions  relevant  to  the  individual  satellites  making  up  each  collaborative 
group  (e.g.,  what  burns  should  be  executed  to  achieve  a  certain  level  of  required  coverage).  This 
decomposition  is  natural  for  satellite  operations. 
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Figure  4:  Three  Tiered  EPOS  Hierarchy 


Another  decomposition,  one  based  around  seven  key  levels  of  decision-making  is  describe  in  Table  3-6.  The 
longer  term,  EPOS-wide  issues  are  determined  at  the  upper  levels,  while  shorter-term  decisions  are  lower 
down  in  the  list.  This  is  a  decision-centric  approach  to  decomposition  compared  with  the  entity-centric 
approach  of  the  three  tiers.  The  planning  and  control  architecture  modules  of  Figure  2  can  be  applied  within 
each  of  these  decision  levels.  Note  that  the  description  of  the  decision  levels  shown  in  Figure  5  presents  the 
capabilities  of  a  complete  EPOS  system,  including  functionality  that  has  not  yet  been  implemented. 
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Figure  5:  Decision  Levels 

The  development  of  EPOS  is  an  evolving  process.  Previous  versions  of  EPOS  [1]  focused  on  an  initial  version 
of  the  observation  level  (which  included  a  person-in-the-loop)  and  an  autonomous  implementation  of  the 
maneuver  level.  The  current  version  of  EPOS,  2.0,  is  focused  on  an  autonomous  implementation  of  the 
observation  level  as  well  as  the  sensor  planning  levels.  The  optimization  algorithms  employed  by  each  of 
those  levels  rely  on  data  that  must  be  derived  through  astrodynamic  modeling.  That  modeling  is  described  in 
the  next  subsection. 


Astrodynamics 

Observation  planning  algorithms  require  knowledge  of  the  evolution  of  target  illumination  and  of  sensor-to- 
target  relative  position.  Dynamic  models  are  used  to  answer  the  following  essential  questions  for  each  time  t 
in  the  planning  interval: 

•  Which  targets  are  in  sunlight? 

•  Which  satellites  can  potentially  view  a  given  target,  in  the  sense  that  the  line  of  sight  from  satellite  to 
target  is  not  obstructed? 

•  How  must  the  sensor  on  such  a  satellite  be  oriented  to  have  the  target  at  the  center  of  its  field  of  view? 

•  Paying  attention  to  pointing  angle  limits,  can  the  sensor  be  oriented  so  that  the  target  is  somewhere  in  its 
field  of  view? 

Write  r(t)  for  the  earth-centered  inertial  (ECI)  coordinates  of  an  earthbound  target  at  time  t,  and  v(t)  for  a  unit 
vector  normal  to  the  earth's  surface  (local  vertical)  at  r(t).  Both  r(t)  and  v(t)  are  determined  from  target 
geodetic  latitude,  longitude,  and  time,  using  an  analytic  representation  for  the  right  ascension  of  Greenwich. 

To  determine  if  the  target  is  in  sunlight  at  time  t,  we  need  ECI  coordinates  u( t)  for  the  direction  to  the  sun.  An 
analytic  representation  is  used.  The  target  is  in  sunlight  when  the  dot  product  v(t )-u(t)  is  positive. 
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Satellite  orbital  motion  is  modeled  using  the  simple  J2  secular  theory,  but  the  software  could  be  easily 
modified  to  use  a  more  sophisticated  orbit  propagator.  ECI  position  coordinates,  s(t),  for  a  satellite  are 
computed  from  its  orbit  elements  using  standard  two-body  mechanics  formulas. 

The  look  direction,  from  satellite  with  coordinates  s(t)  to  target  with  coordinates  r(t),  is  computed  as  follows: 


Lrs(t) 


r(t )  -  s(t ) 
||r(t)-s(f) 


The  satellite  can  potentially  view  the  target,  in  the  sense  that  the  line  of  sight  is  not  obstructed  by  the  earth,  if 
the  following  dot  product  condition  holds: 

v(t)  ■  Lrs  (0  <  0 


Satellite  attitude  must  be  modeled  in  order  to  answer  questions  about  sensor  orientation.  Let  As  (t)  denote  the 
nominal  attitude  of  the  satellite  with  coordinates  s(t)  at  time  t.  The  development  version  of  our  system 
assumes  sun-nadir  steering  for  each  satellite;  one  spacecraft  axis  points  to  nadir  and  a  second  axis  is  normal  to 
the  plane  spanned  by  nadir  and  the  sun  direction.  An  operational  planning  tool  would  support  considerably 
more  flexible  specification  of  satellite  attitude  profiles. 

Attitude  As  (t)  is  understood  to  be  a  coordinate  transformation  from  ECI  frame  to  nominal  spacecraft  body 
frame.  Sensor  orientation  is  specified  as  a  second  coordinate  transformation,  As  (t) ,  from  nominal  spacecraft 
body  frame  to  a  frame  fixed  in  the  sensor.  The  composite  transformation,  As  (t)As  (t) ,  could  be  realized  in 
either  of  two  ways: 

•  Steer  the  spacecraft  so  that  As  (t)  is  the  transformation  from  ECI  to  actual  spacecraft  body  frame.  Gimbal 
the  sensor  so  that  As  (t)  is  the  transformation  from  spacecraft  body  frame  to  sensor  frame. 

•  Steer  the  spacecraft  only,  so  that  As  ( t)As  (t)  is  the  transformation  from  ECI  to  the  sensor  frame. 

In  either  case,  we  regard  As  (t)  as  given  and  As  (t)  as  something  to  be  chosen. 

Let  i  be  a  unit  vector  in  the  sensor  look  direction,  expressed  in  the  sensor  frame.  The  target  is  at  the  center  of 
the  sensor  field  of  view  if  the  following  condition  holds: 

As(t)AJt)LrJt)  =  £ 


This  condition  determines  As  (t)  up  to  rotation  about  £  .  It  determines  the  first  two  angles  x(t)  and  y(t)  in  a 
suitable  Euler  angle  sequence,  in  other  words. 

The  target  remains  somewhere  in  the  sensor  field  of  view  provided  that  these  Euler  angles  are  chosen  from 
intervals  \x(t)~  Ax,  x(t)  + Ax\  and  [y(t)-Ay,  y(t)  +  Ay] ,  respectively.  Here  At  and  Ay  are  constants  of 
the  sensor  (field  of  view  half-widths). 

Four  additional  parameters  of  the  sensor,  xmin ,  xmax,  ymjn  and  ymax  are  used  to  specify  feasible  intervals  for 
the  two  Euler  angles.  For  the  sensor  to  be  able  to  view  the  target  at  time  t,  it  is  necessary  that  the  line  of  sight 
be  unobstructed  (dot  product  condition  given  earlier)  and  that  the  intervals  I x  and  I  be  nonempty,  where 

I ,  =  U(t)  -  Ax ,  x(t)  +  Ax]nk,n ,  xmax  ] 

Iy=[y(t)~Ay,  y(0  +  Ay]n[ymin,  ymax] 
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Observation  Planning 

The  observation  planning  level’s  problem  is  to  find  the  primary  target  for  each  sensor  so  that  the  total  science 
value  at  time  t  (including  extra  for  simultaneous  viewing)  is  maximized  provided  that: 

•  no  more  than  one  of  each  type  of  sensor  is  tasked  for  each  target1, 

•  no  more  than  two  (different)  sensors  are  assigned  per  target 

•  demand  for  data  storage  on  a  satellite  does  not  exceed  available  capacity. 

The  problem  is  specified  through  static  data  that  describes  the  available  satellites  and  their  associated  ground 
stations  and  dynamic  data  that  specifies  that  targets  of  interest,  which  targets  can  be  viewed  by  a  particular 
satellite  and  which  ground  stations  can  communicate  with  a  satellite  at  a  given  time.  In  the  discussion  that 
follows,  the  terms  satellite  and  sensor  are  often  used  interchangeably.  The  schedule  is  being  done  for  each 
sensor.  Conceptually,  one  can  think  of  a  satellite  containing  a  single  sensor.  It  is  of  no  consequence  to  the 
formulation  that  follows  if  multiple  sensors  are  contained  on  a  single  satellite.  To  being  the  derivation  of  the 


problem,  the  following  sets  must  be  defined: 

N 

number  of  satellites/sensors 

M 

number  of  targets 

L 

number  of  ground  stations 

Seo  = 

Set  of  sensors  of  type  EO 

SIR  = 

Set  of  sensors  of  type  IR 

S  SAR  = 

Set  of  sensors  of  typeSAR 

Dynamic  Inputs 

The  dynamic  inputs  to  the  observation  level  optimization  specify  data  that  is  changed  from  one  pass  to  the 
next.  The  set  of  targets  that  a  particular  sensor  can  view  at  time  t  are  specified  as: 


jlif  satellite/sensor  i  can  observe  target  j 
[O  otherwise 


The  set  of  ground  stations  that  a  satellite  can  communicate  with  at  time  t  are  specified  as 

[l  if  satellite/sensor  i  can  communicate  with  ground  station  g 
[0  otherwise 


The  values  of  o]-  and  q’u  are  derived  using  the  methods  described  in  the  subsection  describing  astrodynamics. 

Whenever  one  or  more  sensors  view  a  target,  scientific  value  is  gained  by  the  observation.  This  value  is 
specified  quantitatively  as: 

Ujj  =  value  from  satellite/sensor  i  observing  target  j 

w,-  •  •  =  additional  value  from  satellites/sensors  /',  and  i2  observing  target  j  simultaneously 


This  is  a  simple  model  to  handle  the  assumption  that  two  observations  by  the  same  type  sensor  during  any  time  period 
does  not  have  twice  the  value  of  one  observation. 
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Decision  Variables 

In  time  period  t,  the  decision  variables  are: 


[lif  satellite/sensor  i  observes  target  j 
[O  otherwise 

[lif  satellites/sensors andz2  observe  target  j  simultaneously 
0  otherwise 


Static  Data 

Static  data  is  used  to  specify  the  data  storage  capabilities  of  the  satellites.  This  data  is  static  in  the  sense  that  it 
is  particular  to  a  sensor/satellite  only  and  does  not  depend  on  any  targets  or  their  visibility  status.  There  is  an 
assumption  that  ground  resources  are  sufficient  for  any  number  of  the  satellites  to  be  downloading 
simultaneously.  Thus  each  satellite  picks  the  ground  station  with  which  it  has  the  maximum  transfer  rate. 

K :  =  total  capacity  of  satellite  i  (number  of  observation  units  that  can  be  stored) 

k\  =  capacity  on  satellite  i  at  the  start  of  time  period? 

aE0  =  data  units  per  observation  consumed  by  an  EO  sensor 

oIR  =  data  units  per  observation  consumed  by  an  IR  sensor 

aSAR  =  data  units  per  observation  consumed  by  an  SAR  sensor 

d'  =  maximum  number  of  data  units  during  t  ime  period  t  that  can  be  downloaded 

from  satellite  i  when  a  ground  station  is  visible 

The  values  of  k\  can  then  be  defined  recursively  as: 


Formulation 

Given  the  data  definitions  specified  above,  the  mathematical  programming  formulation  is  as  follows  for  each 
time  period  t. 
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max 


X 


Y,x'ui‘,,o'l  +  2>>,jo'.  +  Y.  XtijUijOtij  + 

?G  ^EO  iG  $IR  iG  SSAR 

Z  y‘^jwhhJo[o'2  +  Z 

V*1G  ^£0  ’*2G *1G  ^£0  ’heSSAR 


i2jW:  ,  .O-  O-  + 
1  J  h>l2’J  h  l2 


Z/M: 

?igZa/?  ’h^SjR 


2  ,jW:  :  i0\0\ 
1  J  h’h  J  h  l2 


Subject  to 


Z  x'j 

for  i  =  1  to  N 

j=\..M 

z^ 

<i 

for  i  =  lto  N 

7=1  to  M 

Z  x'ij 

<i 

for  /  =  1  to  AT 

iG  ^EO 

Z  x*v 

<i 

for  j  =  1  to  M 

ie  SIR 

<i 

for  /  =  1  to  AT 

IG  S  SAR 

Z  4 

<2 

for  /  =  1  to  AT 

1=1  to  N 

hi  +  x  hi  ~2y  it,i2j 

>0 

for  tj  A  i2 ,  /  =  1  to  Af 

hi  +  x  hi  “2 y  i, ,i2,j 

<1 

for  tj  *  i2 , 7  =  1  to  M 

t 

X  ij 

e  {0,1} 

for  i  =  1  to  N,  7  =  1  to  M 

I 

y  h’h’j 

e  {0,1} 

for  i  =  1  to  N,  7  =  1  to  M 

The  output  of  the  top-level  optimization  is  at  each  time  f  and  for  each  sensor  the  primary  target  it  is  tasked  to 
look  at  and  whether  the  satellite  downloads  observation  data  or  not. 


Sensor  Planning 

The  assignments  of  targets  to  sensors  at  each  time  period  are  input  to  the  sensor  optimization  level.  The  output 
is  a  pointing  command  for  each  sensor  for  each  time  period.  The  sensors  are  modeled  as  being  mounted  on 
platform  with  two  orthogonal  gimbals.  These  gimbals  are  referred  to  as  the  x  and  y  gimbal  angles. 

Data 

To  begin  the  definition  of  the  data  required  for  this  optimizer,  first  define  the  following  sets: 

R  =  the  set  of  targets 
T  =  the  time  horizon  for  the  problem 

The  model  of  the  sensor  includes  constraints  on  the  movement  of  each  of  the  gimbals.  These  constraints 
include  both  bounds  on  the  gimbal  angles  and  slew  rate  constraints.  The  following  data  is  defined  for  each 
t  =  \..T. 


Bit ) 

Ar  (t).[lu][xy] 

P.[lu]  [xy] 
5.[xy] 


{re  R  |  r  is  assigned  to  the  sensor  at  time  t} 

required  bounds  on  the  gimbal  angles  for  the  sensor  to  be  able  to  see 
target  r  e  Bit)  at  time  l  (e.g.  Ar  ( t).ux  is  upper  bound  on  x  direction) 
gimbal  angle  limits  for  the  sensor  in  each  direction  (rad) 
slew  rate  limits  for  the  sensor  in  each  direction  (rad/At) 


The  notation  [xy]  is  used  to  represent  combinations  of  the  specified  variables.  For  example,  S.[xy]  indicates 
that  both  variables  5.x  and  S.y  are  defined. 
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Decision  Variables 

The  primary  decision  variables  are: 

p(t).[xy]  = 

w,.  ( ? ) 


girnbal  angle  for  the  sensor  at  time  ?  for 
each  direction 

f  1  if  the  sensor  can  see  target  r  e  B(t)  at  time  t 
0  otherwise 


The  variables  are  defined  for  all  t  —  I..T  and  the  values  of  p(0).[xy]  are  specified  as  input.  The  wr(t)  variables 
are  obviously  a  function  of  the  p{t)  variables.  Specifically,  wr  (?)  =  1  if  and  only  if  the  following  are  true: 

A,.  ( t).lx  <  p(t).x  <  A,.  ( t).ux 
A,,  (t).ly  <  p(t).y  <  A,  (t).uy 


In  support  of  this,  define  the  following  binary  variables: 

lif  Ar(t).lx  <  p(t).x 


br(t).lx  = 
br(t).ly  = 

The  relationship  between  wr  (?)  and  br  (?)  is  then 


0  otherwise 

lif  Ar(t).ly  <  p(t).y 

0  otherwise 


br(t).ux  = 
b,(t)-uy  = 


flif  p(t).x  <  Ar(t).ux 

1 0  otherwise 

flif  p(t).y  <  Ar(t).uy 

[0  otherwise 


4wvr  (?)  <  bvr  (t).lx  +  bvr  (t).ux  +  bvr  (t).l.y  +  bvr  ( t).uy 

The  br  (?)  valuables  are  set  with  the  following  constraints 

p(t).x  >  Ar  ( t).lx-M(l-br  ( t).lx ) 
p(t).x  <  Ar  {t\ux  +  M  (1  —  b r  ( t).ux ) 
p(t).y>Ar(t).ly-M(l-br(t).ly) 

P(t).y  <  Ar  ( t).uy  +  M(\-br(t).uy ) 

where  M  is  a  larger  than  the  absolute  value  of  any  A,.  (?)  entry. 


Constraints 


The  pointing  limits  for  the  sensor  are  specified  simply  as: 

P.lx  <  p(t).x  <  P.ux  V?  =  1..T 

P.Iy  <  p(t).y  <  P.uy  V?  =  1..T 


The  slew  rates  arc  specified  as: 


| p(t).x- p(t-l).x\  <  S.x 
p(t).y-p{t-\).y\<S.y 


V?  =  \.T 
V?  =  \.T 


Objective  Function 

The  objective  function  maximizes  the  value  gained  from  achieving  all  of  the  assignments.  This  is  calculated  by 
multiplying  the  wr(?)  variables  by  Vr(t) . 
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The  null  position  of  the  sensor  is  pointing  at  the  long-term  targets,  unless  there  is  value  to  be  gained  from  an 
ephemeral  target2,  the  sensor  should  tty  to  point  in  the  null  position.  To  address  this,  a  cost  will  be  attributed 
to  moving  the  sensor  from  the  null  position.  Specifically,  a  term  will  be  added  to  the  objective  function  to 
minimize  the  absolute  value  of  the  p(t).\xy]  variables.  Define  the  non-negative  variables  Dt+_1(t).[xy]  as: 

D+ (t).x- D~  (t).x  =  p(t).x 
D  +  (t).y-D-(t).y  =  p(t).y 

The  objective  function  will  include  a  term  to  minimize 

T 

y  {p+ (t).x  +  D~(t).x  +  D+ (t).y  +  D~  ( t).y ) 

t= 1 


Formulation 


Given  the  data  definitions  specified  above,  the  mathematical  programming  formulation  is  as  follows. 


max  ^ 

<= i 

Subject  to 


X  Vr(t)wvr  (t)  -  y{p  +  (t).x+  D  (t).x  +  D  +  (t).y  +  D  ( t).y  ) 


re  Bv  (t) 


p(t).x  -  p(t  -  l).x  <  S  .x  V  t  =  1  ..T 

p(t).x  -  p(t  -  \).x  >  -S  .x  V  t  =  1  ..T 

P(t).y  —  p{t  —  \).y  <  S  .y  \/t  =  l..T 

p(t).y  -  p(t  -\).y  >  —S.y  \/t  =  l..T 


p  ( t).x  -  Mb  r  (t).lx  >  Ar  (t).lx  -  M 
p  ( t).x  +  Mb  r  (t).ux  <  Ar  (t).ux  +  M 
p{t).y  -  Mb  r  (t).ly  >  Ar(t).ly  -  M 
P(t).y  +  Mb  r  (t).uy  <  Ar(t).uy  +  M 


Vr  g  Bv(t),t 
Vr  g  Bv  (f),  t 
Vr  e  Bv  (f),  t 
Vr  e  Bv(t),t 


1  ..T 
1  ..T 
1..T 
1  ..T 


4 w r  (f)  -  br  (t).lx  -  br  (t).ux  -  br  (t).ly  -  br  (t).uy  <  0 


Vr  g  Bv  (f),  t  =  1  ..T 


D  +  {t).x-D  (t).x  -  p(t).x  =  0  V t  =  1  ..T 

D  +  {t).y  -  D  (t).y  -  p(t).y  =  0  V  t  =  1  ..T 


P(t).x  g  [P  .lx  ,  P  .UX  ] 
p{t).y  g  [P.ly,P.uy  ] 
D[^](0-[xy]G  91 

wr(f)e  {0,1} 
br(tUlu][xy]  g  {0,1} 


The  value  from  observing  an  ephemeral  target  is  assumed  to  be  greater  than  observing  a  long-term  target. 
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Future  Work 

The  development  of  the  EPOS  system  has  progressed  so  that  both  maneuvering  and  coasting  satellites  can  be 
optimized  to  achieve  optimal  target  viewing.  A  limitation  of  the  current  system  is  that  maneuvering  satellites 
treat  their  sensors  as  having  a  fixed  attitude  while  coasting  satellites  utilize  an  autonomous  sensor  optimization 
level.  The  next  version  of  EPOS,  currently  being  developed,  will  merge  these  two  approaches.  Specifically, 
both  maneuvering  and  coasting  satellites  will  be  addressed  and  all  satellites  will  optimize  the  pointing  of  each 
sensor.  In  effect,  the  current  version  of  EPOS  and  the  capabilities  of  the  previous  version  [1]  will  be  merged 
into  a  single  approach. 
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